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1. GRANT MISSION 


This NASA Grant, 5-1 1575, has covered the active participation of the Principal Investigator, 
Prof. Gordon Pettengill, and his Co-Investigator, Peter Ford, in the Mars Orbiting Laser Altimeter 
(MOLA) Experiment, over a period of two years. This participation has included attending team 
meetings, developing data-reduction software algorithms, and processing data, as well as presenting 
several oral reports at scientific meetings. The research has concentrated on the various types of 
Martian clouds that were detected by the laser altimeter, with results summarized below. 


2. TECHNICAL SUMMARY 

The Mars Orbiting Laser Altimeter (MOLA) instrument [ 1 ,2] carried aboard the Mars Global 
Surveyor (MGS) spacecraft, has observed strong echoes from cloud tops at 1.064 pm on 61% of 
its orbital passes over the winter north pole (235° < L s < 315°) and on 58% of the passes over the 
winter south pole (45° < L s < 135°). The clouds are unlikely to be composed of water ice since the 
vapor pressure of H 2 0 is very low at the Martian nighttime polar temperatures measured by the 
Thermal Emission Spectrometer (TES) (3), and by an analysis of MGS radio occultations [4|. Dust 
clouds can also be ruled out since no correlation is seen between clouds and global dust storms. 
The virtually certain composition for the winter polar clouds is C0 2 ice. 

Cloud Types : We have constructed a database of MOLA clouds, identified by clusters of 
closely grouped echoes lying significantly above the planetary surface. Each set of cloud echoes 
was plotted by altitude and along-track distance (see Figs. 3-5), inspected visually, and assigned to 
a category based on its morphology and relation to surface features. The location of clouds by 
latitude, sub-solar longitude (L s ), and solar incidence angle are shown in Figs. 1 and 2. 
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Fig. 1: Location of MOLA clouds by latitude (ordinate, in degrees) and sub-solar longitude, L s (abscissa, 
in degrees). 
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Fig. 2: Location of MOLA clouds by solar incidence angle (ordinate, in degrees) and sub-solar longitude, 

Ls. (abscissa, in degrees). Values of solar incidence greater than +90° or less than -90° indicate that the sun 
is below the horizon. Note that virtually no clouds, except for those of the ground-hugging type (colored 
magenta), occur when the sun is above the horizon, as would be expected if they are composed of C0 2 ice. 

The colors used to denote the cloud locations in Figs. 1 and 2 indicate the categories that we 
have assigned to them, as exemplified in Figs. 3-5 which show the altitude (km) and distance- 
along-track (km) of typical lidar echoes for the three most common categories. 

Less common categories comprise dome-shaped clouds (colored light blue in Figs. 1 and 2), 
dense isolated formations not apparently associated with surface features, ground-hugging clouds 
(magenta), flat clouds that form at altitudes of less than 1 km in the spring and autumn dawn, and 
occasional crater clouds (brown), which form within large impact craters. 

The data points in Figs. 3-5 have been color coded according to the time dispersion of the lidar 
echoes: red indicates reflection from a cloud top whose density onset varies strongly over no more 
than 3m of range; green and yellow denote intermediate gradients; while blue echoes are the most 
dispersed in range— from clouds whose density discontinuity is spread over 90m or more. Cloud 
echoes that saturated the MOLA receiver are colored magenta, and surface echoes are colored gray. 

Under conditions of Mie scattering, and assuming unit geometric albedo, the number density of 

8—2 3 

C0 2 ice particles required to produce the observed echo strengths is about 2x10 r per m , where 
r is the particle radius in pm, which we infer to lie between O.lpm and 50 pm. If r were less 
than 0.1//m, Rayleigh scattering would predominate, the cross-section would rapidly drop off as r 4 , 
and no lidar echo would be detectable. If r were greater than 50pm, the particles would fall too fast 
to sustain the observed wavefront patterns. 
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Fig. 3: Echoes from a typical cloud formation that is adjudged to contain propagating waves— green in 
Figs. 1 and 2— which are more common in the northern than in the southern polar regions. Ordinate 
represents altitude (km) and abscissa distance along track (km). 



Fig. 4: Example of clouds associated with topographic features— colored blue in Figs. 1 and 2— with the 
appearance of mountain waves or perhaps snow tails. Ordinate represents altitude (km) and abscissa distance 
along track (km). 



Fig. 5: Example of diffuse clouds— colored red in Figs. 1 and 2— whose lidar echoes show very little 
time dispersion. They are seen only in southern polar latitudes. Ordinate represents altitude (km) and 
abscissa distance along track (km). 

N on-reflective clouds: In addition to reflective clouds, MOLA often receives no echo from 
either cloud or surface, thus implying an intervening absorber. These clouds form a distinct 
category. Their distribution in latitude, season, and solar incidence angle is shown in Figs. 6 and 7. 
Since these clouds are not exclusive to the polar night, and are more common during hemispheric 
dust storms, they may be composed largely of dust or diffuse H 2 0 ice, rather than C0 2 ice. 
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Fig. 6: Location of non-reflective MOLA clouds by latitude (ordinate) and sub-solar longitude, Ls 
(abscissa). Those clouds located within Hellas crater are shown in red. 
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Fig. 7: Location of non-reflective MOLA clouds by solar incidence angle (ordinate)and sub-solar 
longitude, Ls. (abscissa). Values of solar incidence greater than +90° or less than -90° indicate that the sun 
is below the horizon. Those clouds located within Hellas crater are shown in red. 

Buoyancy Waves. We have identified and measured all lineaments within the cloud-top echoes, 
such as those visible in Figs. 3 and 4. We assume that the multiple wave fronts of the propagating 
class (Fig. 3) are buoyancy waves moving through a C0 2 atmosphere [5|. Particles of dry ice are 
alternately condensed and evaporated as the ambient temperature is modulated above and below the 
local freezing point by the passing wave. This implies that the atmosphere is close to its “wet” CO 
adiabat at ~148K, at which point the buoyancy frequency is estimated to be ~9.3 1x10 rad/sec. 

Spatial-Temporal Correlation. The combination of the MGS sun-synchronous orbital period, 
~7061 s, and the length of the Martian day, ~88775 s, results in a near overlap of surface tracks at 
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7-day (88 orbit) intervals. However, since the tracks culminate (converge) at ±86° latitude, more 
frequent overlaps occur at polar latitudes. We therefore auto-correlated the locations of all clouds of 
a given category, assuming that they circle with the ambient wind around the pole at a nominal zonal 
flow. Fig. 8 shows the result of these correlations around the south pole, for a range of assumed 
zonal rotation rates. 

Note the strong correlation at non-zero rotation rate for the class of propagating clouds ( e.g ., 
Fig. 3) and the relatively weak correlation between clouds related to topographic features (e.g.. 
Fig. 4), suggesting that the former are moving in a polar vortex — ca. 5 m/sec at 75°S 
latitude— while the latter are not. The tracks of some correlated propagating clouds are shown in 
Fig. 9. 



Fig. 8: Correlation between pairs of clouds detected within 48 hours of each other and separated by > 1° 
of arc measured from the center of Mars, as a function of zonal vortex wind drift rate (degrees of longitude 
per day). Green: correlation of propagating clouds; blue: of topographic clouds; black: of all cloud 
categories. 



Fig. 9: Track history of correlated clouds near the North and South Poles of Mars. Color denotes the 
interval over which correlation was observed: pink: 1 day; green: 2 days; blue: 3 days; cyan: 4 days. 
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Infrared Emission. Simultaneous observations were made by MOLA and by the MGS Thermal 
Emission Spectrometer (TES) [6]. At the low temperatures of the polar night, the emission 
spectrum T x between 15//m and 25// m is most useful. With minimal contribution from H z O and 
dust, the emission is dominated by the C0 2 rotation band at \5pm and by the thermal (Planck) 
spectrum from the surface and atmosphere at their respective temperatures. Although the 
fundamental spatial resolution of each of the 6 TES detectors is ~3x3 km, all six are usually 
summed together for polar nighttime observations, in order to achieve detectable signal levels. Each 
footprint is therefore ~6x9 km and encloses the ~20Q m diameter MOLA footprint. The surface 
emissivity is typically ~0.95 at these wavelengths, so most TES observations in the polar night (see 
Fig. 10) show an atmosphere and surface at -148K. 

Polar “Cold Spots”: Anomalously low emission temperatures have been observed in the polar 
night, however, both by TES and by the IRIS instrument on Mariner 9 [7). Fig. 10 shows TES 
spectra at three latitudes. On the left, at 60°S, the warm sunlit atmosphere emits strongly at \5pm, 
but is almost transparent at other wavelengths, which are dominated by surface emission. At 85.5°S, 
the atmosphere is colder than the surface and absorbs at 15 pm, while the radiometric temperature 
drops at X> 20pm, indicating an anomalous “cold spot.” At 79.5°S, the long-wavelength anomaly 
disappears and only the 15pm absorption feature remains. 

175K 


150K 


125K 


100K 

Fig. 10: Example of a cold spot. The blue line in the image on the left traces the TES track for orbit 
8799, L s = 122°. The figure on the right shows portions of the TES spectra at three locations marked X 
along the orbit track. At 60°S, the upper atmosphere is warm (elevated temperatures in the 15/<m C0 2 
absorption band) and the spectrum beyond 18//m is consistent with 148K physical temperature and 0.95 
surface emissivity. At 85.5°S, a cold spot appears, i.e., the radiative temperature falls below 140K at 
wavelengths longer than 18/<m. At 79.5°S, the anomalously low radiative temperature beyond 18/im has 
disappeared. The red ticks mark the wavelengths at which three contributions dominate: 15/<m: the 
atmosphere above ~10km; \7fim: the atmosphere below ~5km; and 25//m: the surface (physical 
temperature x emissivity). 

Three explanations for the cold spots have been advanced. They could result from C0 2 ice clouds 
[81, from a deep layer of CO, frost [8,9], or from surface layers of C0 2 slab ice [ 10|. In each case, 
the low value of emissivity results from volume scattering, but under very different conditions. To 
choose among these models, we have made maps of the spectral slope, dE/dX,for 20pm < 
X < 25pm at intervals of 5° in L s . The cold spots show up as isolated minima, as shown in Fig. 1 1. 
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Fig. 12: The location of the 163 cold spots identified in TES observations from October 1997 to June 
2001 are shown as white dots in these MOLA topographic maps. In the north, cold spots are found close to 
the eastern rims of 6 of the 20 largest impact craters northward of 55° latitude, and others are clustered on 
east-facing slopes within the permanent polar cap, but there is no significant correlation between cold spots 
and surface altitude. In southern polar regions, there is no significant correlation between cold spots and 
either altitude or surface slope or class of surface feature. As an example, the blue dots in the south polar 
map show the distribution of “spider” formations taken from the work of Piqueux et al. [11] The spot 
locations also seem to be uncorrelated with several other classes of feature that are generally assumed to 
result from the deposition and evaporation of the seasonal cap, e.g., “dalmatian” spots, “fried egg” haloes, 
etc., etc. In the absence of any strong correlation with feature class, we are led to conclude that the 
mechanism that is responsible for the cold spots is restricted to the seasonal CO z layer itself. 

Automatic cold spot identification : We have developed the following process for locating and 
measuring cold spots in TES data and correlating their occurrence with MOLA measurements: 

1. Examine all TES spectra taken within 10° of nadir, and orbits for which MOLA echoes are 
available (3/27/98-6/29/01). Apply the Planck equation to derive radiative temperatures for 


8 




wavelengths, and estimate the average by least squares fitting. Create gridded images of the 
spectral derivatives, ffT^/d'k, for 20/im < X < 25/im, at 5° intervals in solar longitude L s . 

2. Smooth each gridded image with a Gaussian filter of 100 km FWHM. Search for regions 
greater than 50 km diameter in which ffTJd'k < -2.0° K//im. Confirm the identification by 
visual inspection and merge duplicate spots from consecutive images. Then compute the 
cold-spot location and diameter. This procedure results in a total of 163 cold spots, as 
shown in Fig. 12, 84 of them northward of 65.9°N latitude and 79 southward of 67. 1°S. 

3. Examine the TES archive again and determine the time at which each of the cold spots was 
first observed, i.e., when TES footprints falling within a spot located in Step 2, smoothed 
by a 50 km low-pass filter, first fell below ffTJdX < -2.0° K///m. 

Search for elevation change: During polar winters, the average altitude of each polar region has 
been observed to increase by 1-2 meters. If cold spots represent regions of recent C0 2 snowfall, it 
is possible that they will receive a significant fraction of their annual deposit within hours or days. 

1. To investigate possible altitude changes associated with cold spots, we have collected the 
MOLA observations that fall within the boundary of each cold spot during a time period 
extending from 30 days before, to 30 days after its first appearance in TES observations. 

2. Treating each cold spot independently, we begin by examining the orbits in which the 
MOLA footprints were reported. Pairs of orbits are selected and their points of intersection 
computed. Three classes of orbit pairs are distinguished: 

a. both orbits occurred before the cold spot was first detected. 

b. neither orbit occurred before the cold spot was first detected. 

c. one orbit was prior to cold spot detection; the other was during or later. 

3. Eliminating those crossing points that lie outside the chosen cold spot, we use the three 
footprints from each orbit that lie closest to the crossing point to estimate the altitude above 
geoid at the crossing point by quadratic interpolation. This results in two altitudes, h t and hj, 
one from each orbit. We call their difference dh ij = h t - h r 

4. We compute the average <dh> and standard deviation a{Sh} of the Sh tj for each of the three 
classes of crossing points. We also compute the standard deviation o{6h n } of the dh jj for 
each orbit n (i.e., n-i or j) in orbit classes “a” and “b”. In this way, we identify those 
orbits that systematically contribute large values to <Sh>, and whose ephemerides 
determined by range- and doppler-tracking are therefore suspect. After removing these 
orbits from the MOLA data set, we iterate the analysis procedure from Step 2 through Step 
4 until no orbit contributes a o{dh n } of more than 5 meters to < dh>. This resulted in the 
removal of 194 from a total of 5327 orbits. 
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5. Cold spots possessing a wide range of topography are identified by their larger pre-cold 
spot values of <dh> and are discarded from further analysis. 1 

6. Table 1 shows <6h>, averaged over all cold spots for which the pre-spot value of <dh> is 
less than 25 cm. No significant altitude change is found to be associated with cold-spot 
formation 

7. We repeated the analysis, Steps 1 through 5, varying the period of pre- and post-cold spot 
formation from 7 to 60 days, and using o{6h} and o{h } to filter the data as described in 
Step 5, above. Finally, we reran the analysis assigning random longitudes to the cold spots. 
In each case, no statistically significant altitude change was detected. 


Table 1. Altitude variance <8h> averaged over all cold spots for which the 
pre-spot value of <8h> is less than 25 cm. 

North South 


Number of Spots 


a. Pre/Pre (cm) 

b. Post/Post (cm) 

c. Pre/Post (cm) 


58 28 

0.1±3.1 -0.8±7. 1 

3.3±5. 1 24.9±10.2 

2.9±4. 1 21.2±18.4 



Fig. 13: Orbit Crossing Points: (a) Orbits prior to first cold spot appearance; (b) Orbits during and after 
cold-spot appearance; (c) One orbit before, the other after first cold-spot appearance; (d) Consecutive MOLA 
footprints used to interpolate the altitude at an orbit crossing-point. Orbit tracks prior to the first appearance 
of the cold spot are shown in red; the remainder in blue. Crossing-points that lie within the cold spot (large 
black circles) are shown as filled black dots; those lying outside are open black circles, and will not included 
in subsequent calculations. The altitude of each valid crossing-point is estimated independently from each of 
the two orbits by interpolating from the nearest triplet of MOLA footprints, (see Figure d, above, in which 
only the post-cold spot triplet, the red circles, is shown.) 


Daylight clouds: While most MOLA clouds were seen during the polar night, and were most 
likely composed of C0 2 ice crystals, some 3% occurred in daylight, mostly at high latitudes. In the 
north, these daytime clouds were restricted to three seasons: late winter (260° < L s < 10°), mid- 
spring (L s ~ 60°), and late fall (140° < L s < 190°). This seasonal pattern was similar to that of the 

1 We could instead have used o{dh} or even o{h}, the standard deviation of pre-spot MOLA altitudes, as a similar 
measure of kilometer-scale roughness. In practice, the three metrics are highly correlated, but <8h> is the least 
sensitive to regional slopes. 
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non-reflective clouds (see above). In the south, daytime clouds were seen in the fall (30° < L s < 70°) 
and late winter (100° < L s < 190°), but at significantly higher latitudes than non-reflective clouds. 
The TES instrument measured thermal emission profiles for 58% of the MOLA daylight clouds. In 
no instance did the emission temperature drop below 150K in the 15—18// wavelength range, 
indicating that these clouds were composed of H 2 0 ice or dust, but not CO z ice. 

Correlation with MOC Images: A distinctive class of daytime cloud was restricted to the late 
northern spring. Its MOLA echoes were closely grouped at low altitude (below 2 km, vs. 5-15 km 
for polar CO z clouds), and showed no evidence of trapped or propagating buoyancy waves. 
Because of the MGS sun-synchronous orbit, these clouds were observed in only two ranges of 
solar incidence angle: 50-55° and > 75°. In each instance, TES showed that these clouds lay at the 
boundary between an isothermal, stably stratified polar atmosphere and an adiabatic atmosphere 
overlying a warmer surface. Whenever they were imaged by the MOC wide-angle camera, the 
clouds were seen to lie at the extreme edge of the seasonal cap, as shown in Fig. 14, and were most 
probably composed of trapped H 2 0 ice and/or dust that was being released as the winter C0 2 
deposits evaporated. 



Fig. 14: Examples of daylight “surficial” MOLA clouds (green) overlayed on MOC images. In orbit 
6500 (left), the cloud overlies the southern extremity of the seasonal cap, which is evaporating in the dawn 
sunlight. The same may be true of orbit 1387 (right), but the edge of the seasonal cap is obscured by a dust 
storm. 

• 3.3% of MOLA reflective clouds were seen in daylight. Of these, 20% were of the rare 
“surficial” type 

o Of these, 55% were within MOC frames 

o Of these, 78% were located within optically dense clouds at the edge of the polar 
hood 

• 5 1% of MOLA “non-reflective” clouds seen were in daylight 

o Of these, 20% were within MOC frames 

o Of these, 53% were located within optically dense clouds at the edge of the polar 
hood or within topographic depressions 
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• All clouds seen by both MOLA and MOC were observed shortly after dawn, never later. 

o In subsequent MOC frames, these clouds often became pure dust storms 

o The preference of the reflective clouds for low altitudes, and of all MOLA daytime 
clouds for the dawn hours, suggests that they are composed of H 2 0 ice condensed 
onto dust grains. 

o This is corroborated by the negative correlation between MOLA surficial cloud 
reflectivity and TES 18// radiative temperature 
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Fig. 15: TES spectra of surficial clouds observed in orbit 6883. The atmosphere at the southern (left) 
edge of the plot has already been warmed by sunlight for ca. 60 minutes, whereas that at the northern (right) 
edge is just emerging into daylight. The effect of the warming is seen in the increasing 17/*m and 25fim 
emission temperatures, while the upper atmosphere observed in the 15/zm CO, absorption band remains 
cold. 
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3. MEETINGS and PUBLICATIONS 

Team Meetings-. 

6 meetings of the Laser Altimeter Team have been held during the reporting period; of these, the 
PI and/or the co-I have attended 5. 

Presentations: 

Mars Atmosphere Modeling and Observations Workshop, Granada, Spain, January 
13-15, 2003 (Ford), paper presented: “Martian Polar Clouds.” 

Columbia Univ., Astronomy colloquium, Mar 26, 2003, (Ford); paper presented: 
"Martian Polar Weather. ” 

Div Plan. Sci. of the Am. Astron. Soc., Monterey, CA, September 2-4, 2003; (Pettengill 
& Ford); paper presented: “Martian Daytime Clouds Observed by MOLA, TES, and 
MOC. ” 

Am. Geophys. Un., San Francisco, December 8-12 2003; (Ford), paper presented:: 
"Change of Martian surface height associated with polar cold spots." 

Other travel: 

On January 28-31 2002, Peter Ford visited Francis Forget at the Laboratoire de Meteorologie 
Dynamique (LMD) du CNRS, Paris to gain “hands on” experience using the LMD computers 
that run the Paris-Oxford general circulation model of the Martian atmosphere. 

Publications. 

Ford, P.G. and Pettengill, G.H., "Martian Polar Clouds", in preparation. 


4. INVENTIONS 


No patentable inventions have resulted from the research supported under this grant. 
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Figure 5 Jovian auroral X-ray spectrums between 300 eV and 1 keV of the north auroral 
region (a, upper panel) and the south auroral region (b, middle panel) as defined in figure 
4. The vertical dotted line at 0.3 keV shows the low-energy cutoff for the Jovian spectra. 
Each spectral point represents 1 0 measured events which are grouped to allow further 
fittings. The lower global count as well as number of events of the south auroral zone 
compared to the north is clearly visible. Lower panel: (c) Chandra ACIS-S spectrums of 
comets Linear S4 [S4] and McNaught-Hartley [MH]; note that spectrum for comet MH is ; 
plotted after scaling by 0.5. Two noticeable features present in Jovian spectra but absent 
in comets are at around 0.65 keV and 0.75 keV. 






